Quantum control of donor electrons at the Si-Si02 interface 
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Prospects for the quantum control of electrons in the silicon quantum computer architecture 
are considered theoretically. In particular, we investigate the feasibility of shuttling donor-bound 
electrons between the impurity in the bulk and the Si-Si02 interface by tuning an external electric 
field. We calculate the shuttling time to range from sub-picoseconds to nanoseconds depending on 
the distance (~ 10 — 50 nm) of the donor from the interface. Our results establish that quantum 
control in such nanostructure architectures could, in principle, be achieved. 
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Silicon based structures are among the most promis- 
ing candidates for the development of a quantum com- 
puter H 0, H S H S due to the existing high level 
of nanofabrication control as well as the well-established 
scalability advantages of Si microelectronics. Differ- 
ent architectures have been proposed in which nuclear 
spin 0,13, electron spin 0] or electron charge are used 
as qubits. All of these proposals ultimately rely on the 
quantum control of single electrons bound to donors. The 
ability to move or 'shuttle' electrons between a donor and 
the Si surface using external electric field is an essential 
element of Si quantum computer architectures because 
the measurement of the electron spin states can only 
occur at surfaces whereas the qubit entanglement takes 
place at the donor sites. In particular, this shuttling time 
must be much shorter than the spin dephasing time (~ 1 
ms in bulk Si). Application of electrostatic potentials 
at surface electrodes would drag the electron from and 
to the donor allowing the manipulation of the electron- 
donor coupling. Quantum control of donor states is also 
a crucial consideration in Si device 'roadmap' as minia- 
turization of transistors leads to only a few dopants per 
device. 

In this Letter we theoretically consider the problem of 
quantum control of donor-bound electrons near a Si-Si02 
interface. We investigate the precise extent to which a 
donor bound electron in the bulk (a few tens of nanome- 
ters from the interface) can be manipulated between the 
donor-bound state and a surface-bound state (within a 
few nanometers of the interface) by suitably tuning an 
external electric field applied perpendicular to the inter- 
face. We address two issues of paramount importance 
in this context: (i) How fast can this electron shuttling 
between the donor and the interface be done in realistic 
Si structures? (ii) Is the shuttled electron at the surface 
still localized in all three dimensions (which will allow 
to take it back to the donor) or is it a delocalized two- 
dimensional electron (which will make it impossible to 



measure its spin)? Our quantitative answers to these 
questions indicate that quantum control consistent with 
Si quantum computer architectures could, in principle, 
be achieved. 

We consider a single electron bound to a substitutional 
P donor in Si, a distance d from an ideally flat Si/Si02 
(001) interface, under an applied uniform electric field 
F perpendicular to the interface. Although conceptually 
simple, this one-electron problem has no formal solution. 
The approach adopted here is based on well established 
approximations pi 1^. Hfl Ull]. proposed in the context of 
conventional Si-based devices, and validated by exten- 
sive studies available in the literature 01 • The formal- 
ism, briefly outlined below, allows us to keep the essen- 
tial physical aspects of the problem within a clear and 
realistic description. The Hamiltonian is written in the 
single- valley effective-mass approximation |l l| : 



H = T + neFz - 



where T 



2Q 



r y/ P 2 + (z + 2d) 2 2(z + d) 



(1) 



giTj^T j, 7 = m x /m\\ is the ra- 
tio between the transverse (m± = 0.191m) and longitu- 
dinal (mil = 0.916 m) effective masses, accounting for the 
Si conduction band valley's anisotropy, and p = (x,y). 
This equation refers to one of the band minima a long 
z, which become lower in energy for this geometry |12| . 
Lengths and energies are given in rescaled atomic units: 
a* = ti 2 ei/m±e 2 = 3.157 nm and Ry* — m±e i /2fi 2 ef = 
19.98 meV respectively, k = 3.89 x 10~ 7 ef (m/m ± ) 2 
cm/kV, the electric field F is given in kV/cm, and 
Q = (e 2 — £i)/(e2 + ei), where t\ = 11.4 and e 2 = 3.8 
are the Si and Si0 2 static dielectric constants. The sec- 
ond term in Eq. (|TJ is the electric field linear potential, 
the third term is the donor attractive potential, and the 
last two terms are the donor and electron image poten- 
tials respectively. In this case Q < 0, meaning that the 
images keep the same sign as their originating charges. 



The overall potential profile for the donor electron 
along the z-axis is schematically shown on the inset in 
Fig. n Note that this is equivalent to an asymmet- 
ric double-well configuration: The well near the inter- 
face is denoted by A and the one around the donor site 
by B. Assuming the wells are not coupled, we obtain 
variationally each well's ground state, defining a basis 
set {^a, ^b} for the low-lying energy eigenstates of the 
double-well. Truncating the Hilbcrt space into this par- 
ticular two-dimensional subspace, based on the double- 
well analogy, is meaningful only for sufficiently large 
donor- interface separations. If d < a*, a single well de- 
scription is more appropriate , and we therefore limit 
the range of distances examined here to d > 2a* . More- 
over, each well's ground state must have an excitation 
gap to the first excited state that is much larger than 
fcflT, which sets an upper bound for the temperature 
as well as for d. For d — > oo, the interface does not 
play a role and the single-valley approximation breaks 
down 13] . Uncoupled effective mass Hamiltonians for 
the donor electron in the A- and S-regions are written 
as Hi = T + Vi, i = A,B, with 



Va = neFz - 



2Q 



^Td 2 ^Td 2 2(* + d)' 



(2) 



and Vb = — 2/r. The donor-related terms [third and 
fourth terms in Eq. |JT]|] are approximated in Va by their 
value at the interface: Vp+ « 2(— 1 + Q)j \J p 2 + d 2 , pro- 
viding confinement along x-y in the A-region. Further 
assuming that d ^> p leads to the 2-D parabolic po- 



i-iJ pa 

tential approximation suggested in Ref. |l4f V^, a rab (p) = 
(1 - Q)(—2/d + p 2 /d 3 ). The barrier at the oxide inter- 
face is assumed to be infinite, so ^ = for z < —d. 
The following properly normalized variational forms were 
adopted for z > -d [3.Hl|: 



f a (z) x gp(p) 



21+1 

a 2 



{z + <£fe- a W* x Ae-W (3) 

(4) 



* B oc {z + d)e-Vp 2 / a2 + z2 / b2 



where a, (3, a and b are variational parameters chosen to 
minimize Ei = (^i\Hi\^i) for i = A, B. In f a (z) we have 
used I = 2 which gives better agreement with the exact 
wave-function for the infinite triangular well than 
I = 1 or £ = 3, as well as a lower variational energy than 
£ = 1, and essentially the same as £ = 3. For d > 2a* , we 
find that a and b coincide with the Kohn-Luttinger (KL) 
variational Bohr radii for the isolated impurity (d — > oo), 
where a = 2.365 nm and 6 = 1.36 nm. 

A relevant question in the present context regards the 
electron confinement parallel to the xy plane when it is 
drawn towards the interface (A-region) by the field. It 
still remains bound to the donor through Vp+ . The calcu- 
lated energy as a function of d, E p = (g\HA\g), is plotted 
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FIG. 1: (Color online) The broad lines, solid and dashed re- 
spectively, give the ground and first excited state energy terms 
(E p and E' p ) obtained variationally as a function of the donor- 
to-interface distance. The narrow lines give the same energies 
calculated within the parabolic approximation. The inset dis- 
plays an outline of the double-well potential. 



in Fig. ^ The energy of the first excited state calculated 
variationally assuming the functional form g' oc x gp> (p) , 
E' p = (g'\HA\g r ), is also shown. For d = 30 nm, we get 
significant binding (\E p \ ~ 5 meV) for the ground state. 
Results obtained within the parabolic approximation are 
also shown in Fig. ^ As expected, this approximation 
underestimates the binding energies and overestimates 
the gap between successive levels, though convergence 
towards the variational results is obtained as d increases. 
The variational parameter /3, characterizing the radial 
confinement of the ground state parallel to the interface, 
is given in Fig. Efa). For the coherent manipulation of 
electrons in quantum devices, it is crucial that the entan- 
glement of the electronic states occurs in a completely 
controlled and reversible manner. This requires that the 
ionized state near the barrier remains laterally bound 
to its respective donor site, setting an upper bound for 
the operating temperature [ksT <C min(\E p \,E' p — E p )] 
as well as for the donor planar density [n < (/3/2) 2 ] to 
avoid significant wavefunction overlap among electrons 
bound to neighboring donors. For d = 30 nm, we get 
an excitation gap E' — E p ~ 1 meV, and donor electron 
wavefunction confinement within a ~ 40 nm diameter 
region parallel to the interface. These parameters yield 
an upper bound of n ~ 10 10 cm -2 for planar donor den- 
sities and limit the operating temperature to a few K. 
Above these limits, the electrons at the Si surface would 
either form a delocalized impurity band (n > 10 10 cm~ 2 ) 
or become thermally excited (T > 1 — 5K). 

The double- well problem is solved through direct diag- 
onalization ofEq. Q, H = T + V A + V B + 2/ ^ p 2 + d 2 , in 
the non-orthogonal basis {^a,^b}, leading to the two 
lowest energy eigenstates and and eigenvalues Eq 
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FIG. 2: (Color online) (a) Confinement lengths l/f3 and l/a c 
(for a calculated at F = F c ) obtained variationally. (b) Ex- 
pectation value of the electron z-coordinate versus electric 
field intensity F for three values of the donor distance to the 
barrier d. The horizontal lines indicate the barrier position 
in each case. The inset shows the eigenvalues Eo and Ei as a 
function of F for d — 11 nm. 



and E\. The last term in the expression for H is added 
to avoid double counting of the donor potential which 
is partially included in Va through Vp+ . The response 
of the electron to an applied electric field is depicted in 
Fig- Gib), where the expectation value (^ol^l^o) is given 
for three values of d. At very low fields the ground state 
is centered around the donor site ($o ~ ®b), and its 
response to increasing fields is strongly dependent on 
d: For the smaller d values, the electron is smoothly 
drawn from near the P + nucleus toward the barrier as 
F increases, while for the larger d the transition is more 
abrupt, and takes place at lower values of F EEEl. We 
define a critical field F c as the field value at which the 
gap is minimized, (E± — 2?o)min = 9mm, characterizing 
the anticrossing point in a (Eo, E\) vs F diagram [see 
inset in Fig. Gib)]. For stronger fields (F > F c ), 'J'o ap- 
proaches ^ a- In this regime, the variational parameter a 
characterizing the decay of the wavefunction along z [see 
Eq.(GJ)] becomes relevant. Fig.GJa) gives the length l/a c 
vs d for a calculated at F c . We note that for the range 
of distances studied here, l/a c is always smaller than 
a*, justifying the double- well approach discussed above. 
For larger values of F, confinement along the field di- 
rection becomes even stronger, with I/a decreasing by 
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FIG. 3: (Color online) (a) The data points give our estimate 
for the characteristic donor ionization time, defined here as 
r = h/gmin, versus F c . For comparison, we present the in- 
verse ionization rates obtained from the semiclassical approx- 
imation (WKB) for the isotropic hydrogenic model with ef- 
fective Bohr radii equal to each of the variational KL param- 
eters a and b. (b) Critical field versus distance of the donor 
to the interface. The solid line gives a phenomenological fit 
F c oc 1/| 20 1- The dashed line is a fitting to tight-binding re- 
sults in Ref. [r?| of the form F c oc 1/d. (c) Donor ionization 
tunneling (solid line) and adiabatic passage (dashed) times 
versus d. 



about a factor of 3 as F increases from F c to 10F C . The 
field-independent confinement length parallel to the in- 
terface, (1//?), is typically one order of magnitude larger 
than l/a c , and both increase sublinearly with increasing 
d. We now focus on a key parameter for general device 
applications: The tunneling time for donor ionization un- 
der an applied field. This question has traditionally been 
addressed in the literature |lj| in analogy with the hydro- 
gen atom ionization problem, based on the semiclassical 
WKB approximation. The expression for the tunneling 
time of an isotropic hydrogenic atom with effective Bohr 



4 



radius a e g in Si under a uniform electric field F is |19| 
« B (F) = g»^exp(3^) . (5) 

We present a fully quantum-mechanical estimate for the 
tunneling time, which we relate here to the anticrossing 
energy gap through the uncertainty relation. Fig. OJa) 
presents our results for r = h/g m m vs critical field, 
allowing direct comparison with the WKB estimates 
t wkb(-^)- Using for a e s the KL variational Bohr radii a 
and b leads, respectively, to lower and upper bounds for 
the calculated r, within a factor of up to 100 in the range 
of electric fields considered here. Our results cannot be 
fitted by the WKB isotropic model with an intermediate 
value of dcff , a result that may be due to the intrinsically 
anisotropic nature of the system, to the final state in the 
ionized regime also being a bound state here, as well as 
to limitations of the semiclassical approximation. 

A summary of our main results as a function of d, 
which is the single fabrication-related parameter in our 
model system, is also presented in Fig.|3J Fig.|3{b) gives 
the critical field versus d results, which are well fitted 
by the function F c oc l/|zo|, where zq — d — 5/ a is the 
saturation value of (z) after complete donor ionization. 
Results varying as F c oc 1/d, obtained through a tight- 
binding model where the Si conduction band details, in 
particular its six- fold degeneracy, are incorporated |17j . 
are also given. The good qualitative agreement between 
the two curves indicates that our approach indeed cap- 
tures the essential physical aspects of the system. The 
quantitative differences may arise not only from the sin- 
gle valley effective mass approximation adopted here, but 
also from the different geometries considered. These dif- 
ferences are manifest in the fact that the relevant dis- 
tances in the phenomenological fittings are |zo| and d 
respectively. Fig. 02c) gives the tunneling times vs d, 
showing that r increases by 5 orders of magnitude as the 
donor distance to the interface increases by a factor of 5, 
with r ranging from sub-picosecond to nanosecond time 
scales (e.g. r ~ 3 ps for d ~ 15 nm). Typical adiabatic 
passage times ^7$ r a — h\e\F c d/ g^^, also shown, are or- 
ders of magnitude larger than the tunneling time (e.g. 
r a ~ 0.4 ns for d ~ 15 nm). 

The results presented here define critical parameters to 
be taken into account in a variety of scenarios where these 
processes are involved. In silicon-based nano-electronic 
devices [i^, the critical field dependence on donor po- 
sitioning and the tunneling time provide relevant infor- 
mation concerning threshold-voltage control and device 
switching times respectively. For the reversible manipu- 
lation of electrons in quantum devices, properties of the 
states bound to the interface define upper bounds for the 
operating temperatures and for the donor planar densi- 
ties required (n < 10 10 cm~ 2 ) to avoid significant wave- 
function overlap among electrons bound to neighboring 
donors. The times required for the shuttling processes 



will be different depending on whether quantum infor- 
mation is stored in spin or in charge degrees of freedom. 
In the case of spin qubits, tunneling times define the 
limiting time scales, since electron tunneling does not 
affect spin coherence. Spin coherence times in bulk Si 
(T 2 ~ 1ms) 21] are at least 5 orders of magnitude longer 
than the tunneling times reported here. Moreover, spin 
coherence in Si can be further enhanced by isotopic pu- 
rification [2ll . Recent experiments |23j | demonstrate 
that T2 is also sensitive to the dopant depth below the 
interface as well as to the interface quality. This means 
that, for the particular geometry of interest here, the Si 
bulk values of T2 give an upper bound for the coherence 
times: Careful interface optimization as well as avoiding 
inhomogeneities in device fabrication and applied fields 
constitute additional requirements for a large number of 
operations to be performed before spin coherence is lost. 
For charge qubits, orbital/charge coherence is required 
and adiabatic evolution of the electron state must take 
place, requiring much longer time scales. Charge relax- 
ation times are much shorter than spin coherence times, 
of the order of 200ns for Si quantum dots surrounded by 
oxide layers |24|, but still longer than the adiabatic pas- 
sage times for d < 20nm. Our realistic results provide 
physical bounds for quantum control of qubits based on 
Si:P donor electron states. In particular, provided that 
spin coherence times near a gated interface are reduced 
by no more than one order of magnitude as compared 
to the bulk values [25|. gate- voltage induced spin ma- 
nipulation on the single electron level may be feasible in 
donor-based Si quantum computer architectures, similar 
to the exciting recent experimental results on single elec- 
tron mani pul ations in gated GaAs quantum dot nanos- 
tructures [26|. 
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